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Abstract — The  bound  electronic  nonlinear  refractive  index,  n2- 
and  two-photon  absorption  (2PA)  coefficient,  6,  are  measured 
in  a  variety  of  inorganic  dielectric  solids  at  the  four  harmon¬ 
ics  of  the  Nd:YAG  laser  using  Z  scan.  The  specific  materials 
studied  are:  barium  fluoride  (BaF2),  calcite  (CaCO. ),  potas¬ 
sium  bromide  (KBr),  lithium  fluoride  (LiF),  magnesium  fluoride 
(MgF2),  sapphire  (AI2O3),  a  tellurite  glass  (75%Te02+  20%ZnO 
+  5%Na20)  and  fused  silica  (Si02).  We  also  report  >>  ,  and  3 
in  three  second-order,  TT  nonlinear  crystals:  potassium  titanyl 
phosphate  (KTi0P04  or  KTP),  lithium  niobate  (LiiNbO;),  and 
/3-barium  borate  (/3-BaB204  or  BBO).  Nonlinear  absorption  or 
refraction  can  alter  the  wavelength  conversion  efficiency  in  these 
materials.  The  results  of  this  study  are  compared  to  a  simple  two- 
parabolic  band  model  originally  developed  to  describe  zincblende 
semiconductors.  This  model  gives  the  bandgap  energy  ( Eg )  scal¬ 
ing  and  spectrum  of  the  change  in  absorption.  The  dispersion  of 
t?  2  as  obtained  from  a  Kramers-Kronig  transformation  of  this 
absorption  change  scales  as  Eg  4.  The  agreement  of  this  theory 
to  data  for  semiconductors  was  excellent.  However,  as  could  be 
expected,  the  agreement  for  these  wide  bandgap  materials  is  not 
as  good,  although  general  trends  such  as  increasing  nonlinearity 
with  decreasing  bandgap  energy  can  be  seen. 

I.  Introduction 

HE  DISPERSION  of  the  nonlinear  refraction  (NLR) 
and  spectrum  of  nonlinear  absorption  (NLA)  are  mea¬ 
sured  in  the  following  samples:  barium  fluoride  (BaFo).  cal¬ 
cite  (CaCOa),  potassium  bromide  (KBr),  lithium  fluoride 
(LiF),  magnesium  fluoride  (MgF2),  sapphire  (AI2O3),  a  tellu¬ 
rite  glass  (75%TeC>2+  20%ZnO+  5%Na20)  and  fused  silica 
(SiOz).1  Measurements  are  made  at  the  four  harmonics  of 
a  TEMoo,  single-pulse,  picosecond  Nd:YAG  laser  (i.e.,  A  = 
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]The  cubic  BaF2  sample  was  supplied  by  Lawrence  Livermore  National 
Laboratory  (LLNL)  and  is  cut  so  that  k  is  parallel  to  [100]  and  E  parallel  to 
[010].  The  tetragonal  MgF2  sample  was  supplied  by  LLNL  and  has  k  parallel 
to  [001]  and  E  parallel  to  [100].  The  cubic  LiF2  sample  was  supplied  by 
LLNL  and  has  unknown  orientation.  The  trigonal  AI2O3  and  trigonal  CaC03 
samples  have  E  parallel  to  a  crystallographic  axis.  The  KBr  sample  was 
obtained  from  Infrared  Optical  Products  of  New  York. 


1064  nm,  532  nm,  355  nm,  and  266  nm).  The  spectral 
behavior  of  nonlinear  effects  in  these  inorganic  solids  is  of 
interest  because  of  their  wide  use  as  windows  and  elements 
in  high-power  laser  systems.  When  these  materials  are  used 
as  passive  elements  in  high-power  laser  systems,  the  intense 
optical  fields  change  the  optical  properties  of  the  material. 
These  intensity-dependent  changes  in  the  material’s  optical 
properties  affect  the  propagation  of  the  incident  light  and 
govern  the  performance  of  the  optical  system.  The  refractive 
part  of  the  nonlinearity  gives  rise  to  self-lensing,  which  can 
lead  to  catostrophic  optical  damage  [l]-[3],  and  self-phase 
modulation  is  becoming  increasingly  important  in  the  design  of 
ultrafast  laser  systems.  The  absorptive  part  of  the  nonlinearity 
can  decrease  the  transmittance  at  high  irradiance  and  can 
provide  a  path  for  optical  damage  in  the  material  at  lower 
irradiances  than  conventional  dielectric  breakdown  [3],  [4], 

In  this  study  of  transparent  dielectrics,  we  are  specifi¬ 
cally  interested  in  the  ultrafast  nonlinearities  associated  with 
two-photon  absorption  (2PA)  and  bound-electronic  nonlinear 
refaction.  The  coefficients  used  for  quantifying  these  processes 
are  the  2PA  coefficient,  0,  and  the  ultrafast  nonlinear  refractive 
index,  no.  The  large  bandgap  energies  of  most  of  these 
materials  allow  measurement  of  no,  and  0  from  the  near- 
infrared  (1064  nm)  to  the  ultraviolet  (266  nm).  These  materials 
are  assumed  to  become  two-photon  absorbing  when  the  photon 
energy',  Ti.uj ,  equals  one-half  the  bandgap  energy,  Eg.  This 
occurs  in  several  of  these  materials  at  the  third  or  fourth 
Nd:YAG  harmonic.  The  data  collected  here  fills  a  relative 
void,  as  there  are  no  previous  studies  where  both  NLA  and 
NLR  were  measured  over  this  large  wavelength  range. 

Over  the  wavelength  range  studied,  we  do  not  observe  the 
large  dispersion  including  sign  change  as  observed  in  smaller 
gap  semiconductors  [5]-[7].  A  simple  two-parabolic  band 
model  was  very  successful  in  predicting  these  nonlinearities 
in  semiconductors  [7],  [8].  It  predicts  the  nonlinear  absorption 
spectrum  and  uses  a  Kramers  -Kronig  analysis  to  obtain  no 
[5] — [1 1].  The  bandgap  scaling  and  frequency  dependence  of 
no  and  0  were  correctly  described  for  semiconductors.  While 
the  agreement  for  semiconductors  was  remarkably  good,  the 
overall  agreement  for  these  wide  bandgap  materials  is  not  as 
good.  Here  the  data  reasonably  follows  the  predicted  scaling 
for  frequencies  below  the  2PA  edge,  however,  for  photon 
energies  near  or  greater  than  the  2PA  edge,  the  data  quickly 
deviate  from  the  predicted  scaling  of  the  two-band  model. 
Moreover,  we  did  not  measure  a  negative  r?,2  for  calcite  at 
266  nm  as  predicted  by  this  model  (also  for  BBO  as  discussed 
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Fig,  1.  Closed  aperture  Z -scans  in  single  crystal  Bald  showing  a 
self-focusing  nonlinearity  at:  a)  1064  nm,  b)  532  nm,  and  c)  355  nm. 
.Solid  lines  are  fits. 

below).  However,  the  bandgaps  used  for  comparison  are 
the  linear  transmission  cutoff  energies,  which  for  dielectrics 
can  be  significantly  below  the  direct  energy  gap  due  to 
indirect  gap  transitions,  impurities,  defects,  and/or  color  center 
absorption.  Using  the  measured  2PA  spectrum  to  yield  an 
effective  direct-bandgap  energy,  E' .  we  find  better  agreement 
with  the  calculated  values  for  r?2  in  the  ultraviolet. 

We  also  report  the  dispersion  of  n-2  and  spectrum  of  2PA 
in  three  second-order  nonlinear  materials:  potassium  titanyl 
phosphate  or  KTP  (KTiOPCL),  lithium  niobate  (LiNb03), 


Fig.  2.  Z-scan  traces  in  BaFo  at  266  nm:  (a)  open  aperture  Z-scan  showing 
2PA,  (b)  closed  aperture  Z- scan  showing  both  nonlinear  refraction  and  2PA, 
and  (c)  divided  Z-scan  showing  a  sclf-focusing  nonlinearity.  Solid  lines  are 
fits. 

and  //-barium  borate  or  BBO  (ii-BaBA), )}  In  these  materials, 
NLA  and  NLR  are  of  importance  since  either  process  can  limit 
the  frequency  conversion  efficiency  (c.g.,  2PA  of  the  harmonic 
lowers  the  efficiency).  Also,  n2  can  alter  the  phase-matching 
conditions  by  changing  the  refractive  index  of  the  material 
in  an  anisotropic  manner.  While  BBO’s  transmission  range 
extends  into  the  ultraviolet,  allowing  measurements  at  up  to 
the  third  Nd:YAG  harmonic,  KTP  and  LiNb03  can  only  be 
measured  up  to  the  second  harmonic  due  to  residual  linear 
losses  at  355  nm. 

II.  Measurements 

We  use  the  Z- scan  [12]  method  that  is  described  in  detail 
in  [13]  to  separately  measure  the  NLA  and  NLR  of  these 

’  The  orthorhombic  KTP  samples  were  provided  by  DuPont.  Two  samples 
were  used,  one  .>■  cut  used  at  532  nm  with  E  parallel  to  a  and  the  other 
z  cut  used  at  1.064  / im  with  E  parallel  to  y.  The  hexagonal  BBO  sample 
was  provided  by  Crystal  Technology  and  had  k  parallel  to  the  c  axis.  The 
hexagonal  LiNb03  sample  was  provided  by  Crystal  Technology  and  is  x  cut 
with  E  parallel  to  the  optic  axis. 
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TABLE  I 

Summary  of  the  Experimentally  Determined  3  and  n2  Along  With  the  Bandgap  Energies  Used  as  Discussed  in  the  Text 


E, 

e’« 

nj  x  10 

44  esu 

P  cm/GW  | 

Material 

(eV) 

(eV) 

1064  nm 

532  nm 

355  nm 

266  nm 

532  nm 

355  nm 

266  nm 

LiF 

11.6 

NM 

2.5±0.5 

1.9±0.4 

1.9±0.4 

4.0±1.0 

=  0 

sO 

=  0 

MgF2 

11.3 

NM 

1.9L0.4 

1.9±0.4 

2.2±0.4 

5.0±1.0 

=  0 

=  0 

£0 

BaF2 

9.2 

9.0 

5.0±1.0 

7.5*1. 5 

9.7±1.9 

11*2 

£0 

£0 

0.06*0.015 

Si02 

7.8 

9.0 

7.4L1.5 

7.8±1.6 

8.5±1.7 

28*6 

sO 

£0 

0.05*0.01 

AhO, 

7.3 

8.7 

13*3 

14*3 

16*3 

26*5 

sO 

=  0 

0.09*0.02 

BBO 

6.2 

R 

OO 

kd 

11±2 

21±4 

14*3 

1±0.3 

sO 

0.010*0.002 

0. 9*0.1 

KBr 

6.0 

NM 

29*6 

47*9 

NM 

NM 

=  0 

NM 

NM 

CaC03 

5.9 

“8.0” 

11*2 

11*2 

14*3 

46*9 

=  0 

0.018*0.004 

0.8±0.2 

LiNbOn 

3.9 

4.4 

48*7 

440*70 

NM 

NM 

0.38*0.08 

NM 

NM 

KTP 

3.8 

4.6 

100*20 

98*15 

NM 

NM 

0.10*0.02 

NM 

NM 

Te  Glass 

3.6 

4.3 

81±16 

430*90 

NM 

NM 

0.62*0.12 

NM 

NM 

materials.  The  Z- scan  is  particularly  useful  for  this  spectral- 
dependence  study  since  it  is  an  absolutely  calibrated  method 
allowing  measurements  at  different  wavelengths  to  be  com¬ 
pared.  The  sample  thicknesses  and  irradiances  used  allow 
the  “thin”  sample  approximation  to  be  used,  i.e.,  linear  or 
nonlinear  phase  propagation  effects  within  the  sample  are 
small  [14], 

The  Q-switched  and  modelocked  NdiYAG  laser  used 
for  these  experiments  produced  TEMo,o  mode,  single  30- 
ps  (FWHM)  pulses  (switched  out)  that  were  measured  by 
a  second-order  autocorrelation  technique  [15],  To  measure 
the  frequency  dependence  of  n2  and  ,6 ,  we  used  nonlinear 
frequency  conversion  of  the  fundamental  of  Nd:YAG  (1064 
nm)  to  its  harmonics  (532  nm,  355  nm,  and  266  nm)  with 
KD*P  crystals.  Best-form  focusing  lenses  (/  ex  150  mm) 
were  used  for  each  of  the  four  harmonic  beams.  The  beam- 
waists,  wo  (HWl/e2M),  and  pulsewidths,  rp  (FWHM),  of 
the  fundamental  beam  and  its  harmonics  are  respectively:  at 
1064  nm,  35  /j.m  and  30  ps;  at  532  nm,  24  /im  and  22  ps;  at 
355  nm,  19  //,m  and  17  ps;  at  266  nm,  16  /rm  and  15  ps.  All 
irradiances  are  quoted  within  the  sample,  i.c.,  Jo  =  (1  —  R)I, 
where  R  is  the  Fresnel  reflection  coefficient,  and  I  is  the 
external  irradiance.  The  overall  experimental  uncertainty  is 
approximately  ±20%.  The  primary  uncertainties  originate 
from  the  determination  of  the  irradiance  distribution  used 
in  the  experiment,  i.e.,  beam  waist,  pulsewidth,  and  energy 
calibration.  All  the  quoted  energies  are  calibrated  against  a 
Laser  Precision  RJ-7620  energy  meter  using  the  RJP-735 
pyroelectric  detector  head.  Relative  errors  are  reduced  by 
using  energy  windows  and  pulsewidth  windows  with  the  data 
aquisition  system.  For  example,  the  pulsewidth  of  each  laser 
firing  is  monitored  by  measuring  the  ratio  of  the  square  of 
the  energy  in  the  fundamental  to  the  energy  in  the  second 
harmonic  and  only  pulses  where  this  ratio  is  maintained 
within  the  range  of  ±10%  are  accepted  [16].  Narrow  band 
interference  filters  (spike  filters)  are  used  on  the  detectors  for 


TABLE  II 

Linear  Index  of  Refraction  Used  for  the 
Materials  Studied  at  the  Wavelengths  Used 


Material 

1064  nm 

532  nm 

355  nm 

266  nm  1 

LiF 

1.3 

1.3 

1.3 

1.3  | 

MgF2 

1.4 

1.4 

1.4 

1.4 

BaF; 

1.47 

1.49 

1.5 

1.5 

SiC>2 

1.48 

1.48 

1.48 

1.5 

AfA 

1.76 

1.78 

1.8 

1.82 

BBO 

1.6 

1.61 

1.64 

1.66 

KBr 

1.53 

1.55 

CaC03 

1.59 

1.6 

1.6 

1.62 

LiNbO, 

2.2 

2.23 

KTP 

1.77 

1.78 

I  Te  Glass 

2.0 

2.0 

each  Z-scan  measurement  at  a  specific  harmonic.  This  ensures 
the  changes  in  transmittance  monitored  are  not  affected  by 
flashlamp  noise  and  other  spurious  light  sources  that  may  be 
present. 

We  show  example  data  for  BaF2  in  Figs.  1  and  2.  The  linear 
transmittance  of  BaF2  cuts  off  at  ~135  nm.  This  corresponds 
to  an  optical  gap  of  Eg  ~  9.2  eV,  therefore,  BaF2  should 
exhibit  2PA  for  photon  energies  hui  such  that  4.6  eV  <hui  < 
9.2  eV.  This  corresponds  to  wavelengths  <270  nm,  which  can 
be  reached  by  the  fourth-harmonic  of  Nd3+:YAG  at  266  nm. 
Closed  aperture  ±-scans  (linear  transmittance  of  the  aperture 
~40%)  for  BaF2  at  the  fundamental,  second  harmonic,  and 
third  harmonic  of  the  Nd3+:YAG  are  shown  in  Fig.  1.  No 
nonlinear  absorption  was  observed  at  these  wavelengths  using 
open  aperture  Z-scans.  For  each  of  the  Z-scan  traces,  we 
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hco/E 

g 

Fig.  3.  Plot  of  the  sealed  n2  data  as  a  function  of  hu/Eg,  along  with  the  theoretical  dispersion  function  (7 2  (solid  line). 


see  a  positive  n2.  For  example,  using  an  irradiance  I0  ~  103 
GW/cm2  at  532  nm  [Fig.  1(b)],  a  fit  to  the  data  gives  rt2  ~ 
(2.1±0.4)x  1(T16  cmAV  or  (7.5±1.5)x  10-14estt. 

For  an  incident  wavelength  A  =  266  nm,  which  gives  a 
photon  energy  of  Tiu  ~  4.7  eV,  2PA  is  observed  in  BaF2. 
Fig.  2(a)  shows  the  open  aperture  Z-scan  trace  from  which  we 
determine  a  value  of  ft  ~  (0.06±0.015)  cm/GW  for  Iq  ~72 
GW/cm2  The  closed  aperture  Z- scan  is  shown  in  Fig.  2(b). 
The  fit  gives  a  value  of  n2  ~  (3.1  ±  0.6) x:  10“ 16  cm2/W  or 
(1.1  ±  0.2)x  10_13esu.  Dividing  the  data  of  Fig.  2(b)  by  the 
data  of  Fig.  2(a)  gives  data  showing  only  the  refractive  part  of 
the  nonlinearity  as  plotted  in  Fig.  2(c)  [13],  This  figure  shows 
the  signature  of  a  positive  n2.  The  dip  in  the  data  (deviation 
of  the  data  from  the  fit)  arises  from  surface  irregularities  or 
bulk  inhomogeneities  and  is  independent  of  irradiance,  i.e.,  is 
a  linear  background  effect.  This  can  normally  be  subtracted 
by  performing  a  Z- scan  at  very  low  irradiance  levels  where 
the  nonlinear  signal  is  small  [13].  Our  detector  was  not 
sensitive  enough  at  this  wavelength  to  accurately  subtract  this 
background,  and  the  subtraction  was  not  performed. 

The  measured  values  of  n2  and  ft  for  the  materials  studied 
are  summarized  in  Table  I.  Details  of  material  orientations 
and  suppliers  are  given  in  the  footnotes.  The  orientation 
dependence  (i.e.,  anisotropy)  of  the  nonlinear  coefficients  is 
generally  small,  being  less  than  ±25%  and,  thus,  does  not 
significantly  affect  our  conclusions  [8] — [10],  [17], 


In  this  paper  we  define  the  nonlinear  index  by  n  =  n0  + 
n2(i?2),  where  ( )  denotes  a  time  average  and  n0  is  the  linear 
index.  The  relation  between  n2  in  esu  and  in  m2/W  is  given 
by,  n2(esu)  =  (cno/407r)  n2(MKS).  The  values  of  the  linear 
indexes  used  are  given  in  Table  II. 


III.  Modeling 

The  predicted  scaling  relation  for  the  degenerate  2PA  coef¬ 
ficient  as  calculated  from  second-order  perturbation  theory  is 
given  by  [18],  [191 

3(u)  =  R^F2(h^/Eg).  (1) 

n0  rja 


The  function  1ft  depends  on  which  states  are  optically 
coupled  and  gives  the  spectrum  of  degenerate  2PA.  This 
function  contains  all  the  band  structure  related  information  in 
(1),  and  for  the  simplest  possible  case  of  two-parabolic  bands 
is  given  by  [6],  [18],  [19] 


F2(X)  = 


(2x  -  I)3/2 


[2xf 


(2) 


The  Kane  momentum  parameter.  Ep.  has  a  value  near  21 
eV  for  a  number  of  materials  [20],  It  has  not  been  measured 
for  the  materials  we  study  here,  and  we  assume  it  has  this 
same  value  of  21  eV.  Fitting  the  semiconductor  data  of  [19] 
results  in  a  best  fit  for  the  value  of  K  ~  3100,  where  ft  is  in 
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Fig.  4.  Plot  of  scaled  n 2  data  from  previous  measurements  along  with  G 2  (solid  line).  Data  obtained  from  [21]. 


cm/GW  and  Ep  and  Eg  are  in  eV.  We  use  this  value  to  scale 
the  experimental  2PA  data. 

As  shown  in  [5]-[7],  [11],  the  nonlinear  refactive  index  is 
given  by  a  nonlinear  Kxamers-Kronig  integral  of  the  nonde¬ 
generate  nonlinear  absorption  spectrum.  The  scaling  relation 
for  the  degenerate  n 2  is  then  given  by 

K\/Ev 

ri2(esu)  =  - -~~G2{%io/Eg)  (3) 

riQ&g 


where  the  function  G2  contains  the  band-structure-dependent 
information,  and  the  constant  K'  is  directly  related  to  the 
value  of  K.  With  Ep  and  Eg  in  eV,  a  value  of  K  =  3100 
gives  a  value  of  K'  =  7.33  x  10-9  which  is  the  value  used 
throughout  this  paper.  In  [6],  a  value  two  times  larger,  K1  = 
1.5  x  10— 8 ,  was  given.  There,  the  factor  of  two  difference 
between  degenerate  and  nondegenerate  index  changes  was  not 
properly  accounted  for  (see  [9]  for  corrections). 

In  general,  the  Kramers-Kronig  integral  relates  the  nonde¬ 
generate  NLR  to  the  nondegenerate  NLA  and  vice-versa.  Here 
we  give  the  relation  for  the  degenerate  n 2(922)  in  terms  of  an 
integral  over  all  frequencies,  1 o\ ,  of  the  absorption  change  due 
to  a  fixed  input  at  energy  hoj2-  The  dispersion  function  G2  is 
then  given  by 


G 2(22) 


2  r°°  F2(xi;x2) 

*  Jo  xl~  x2 


dx  1 


(4) 


where  Fo(x\ ;  m2)  gives  the  nondegenerate  nonlinear  absorp¬ 
tion  spectral  dependence  and  x\  =  huii/Eq  [6]. 

While  the  functional  dependence  of  the  degenerate  2PA 
coefficient  is  relatively  simple  within  the  two-parabolic  band 
model,  i.e.,  Fo{x),  the  corresponding  dependence  of  the  non¬ 
degenerate  2PA  and,  in  turn,  is  complex  [6],  [9],  [10],  This 
is  further  complicated  by  the  fact  that  NLR  depends  upon 
all  nondegenerate  nonlinear  absorption  processes  including 
electronic  Raman  and  AC-Stark  contributions,  which  must  be 
included  in  ;c2).  The  specific  functional  dependencies 

are  given  in  [6]  for  the  two-parabolic  band  model  (see  [9]  for 
corrections).  There,  the  expected  E~ 4  scaling  in  (3)  is  seen 
as  predicted  by  Wherrett  [18]. 

Fig.  3  shows  our  measured  values  obtained  for  n2  scaled 
by  noEg/ K' sj Ep  versus  Tioj/Eg.  From  (3),  this  gives  an 
experimental  equivalent  of  the  dispersion  function  Go  which 
is  also  included  on  the  same  graph  for  comparison.  Previous 
712  measurements  for  these  materials  are  shown  with  the  same 
scaling  in  Fig.  4  [21] . 

As  seen  in  Fig.  3,  data  agree  with  the  calculation  for  tux 
below  the  2PA  edge  (hut  =  0.5)  within  approximately  a 
factor  of  two.  For  example,  LiF  and  MgF2  possess  the  largest 
band  gaps  (cl  1. 1  eV)  of  the  materials  studied.  Here,  we  see 
relatively  good  agreement  with  the  scaling  predicted  by  the 
two-parabolic  band  model  even  for  the  355-nm  and  266-nm 
data.  MgF2  is  a  tetragonal  crystal  and  has  an  indirect  gap. 
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Fig.  5.  4.  A  log-log  plot  of  the  scaled  n2  data  of  Fig.  3  versus  E,, .  showing  the  predicted  E~A  dependence  (solid  line). 


The  data  for  MgF2  follows  the  scaling  but,  being  somewhat 
smaller,  does  not  agree  in  absolute  magnitude  as  well  as  data 
for  BaF2  or  LiF  which  are  cubic  direct-gap  crystals  [22],  For 
materials  with  smaller  bandgap  energies,  the  experimental  data 
appear  to  diverge  from  the  prediction  of  the  two-parabolic  band 
model  when  the  incident  photon  energy  is  near  the  2PA  edge. 
In  addition,  we  did  not  see  the  predicted  negative  no  values 
for  calcite  or  BBO  at  266  nm. 

Hidden  in  Fig.  3  is  the  several  orders-of-magnitude  variation 
of  n2  due  to  the  i?“4  scaling.  This  variation  is  better  shown  on 
the  log-log  plot  of  Fig.  5.  This  figure  shows  n2  data  multiplied 
by  no  and  divided  by  the  theoretical  dispersion  function  G2 
as  a  function  of  Eg.  Theory  predicts  that  both  graphs  (Figs.  3 
and  5)  are  universal  curves,  independent  of  material.  Thus, 
in  Fig.  5,  data  at  different  wavelengths  for  a  single  material 
should  all  lie  on  top  of  one  another.  Note,  a  datum  for  BBO  at 
266  nm  which  gave  a  very  small  but  positive  n2  is  not  shown 
as  it  lies  well  below  the  bottom  of  the  graph. 

It  is  interesting  to  examine  the  spectrum  of  the  experimental 
2PA  data.  Fig.  6  is  a  plot  of  this  data  scaled  by  n^Eg^EpK 
(i.e.  the  experimental  F2)  as  a  function  of  hu}/Eg.  It  is  appar¬ 
ent  from  Fig.  6(a)  that  the  theory  consistently  overestimates 
ft.  This  would  be  possible  if  the  2PA  edges  of  the  materials 
plotted  were  underestimated.  We  calculate  the  2PA  edge  as 
Eg/ 2  where  Eg  is  determined  from  the  linear  transmittance 
cutoff  frequency,  i.e.,  the  ultraviolet  optical  absorption  edge. 


If  we  assume  that  this  observed  absorption  onset  is  due  to 
impurity  levels  [23],  indirect  gap  transitions,  defects,  or  self- 
trapped  excitons  [24],  the  direct  gap  could  be  significantly 
higher  than  the  linear  absorption  cutoff  energy. 

If  we  allow  the  direct  gap  energy  Eg  to  vary,  so  as  to  make  ft 
fit  (1),  we  obtain  a  new  effective  bandgap  energy,  E'g  —  ftEg. 
For  the  materials  studied  where  2PA  is  observed,  this  results 
in  an  increase  in  I  ft,  (ft  >  1)  except  for  BaF2  which  is  shifted 
~2%  downward,  a  negligible  change.  When  we  make  a  shift 
of  the  x  axis  by  ft  in  Fig.  6(a),  we  must  make  a  shift  of  the  y 
axis  by  <  :! .  We  rescale  data  for  fused  silica  (266  nm),  sapphire 
(266  nm),  BBO  (355  and  266  nm),  calcite  (355  nm  and  266 
nm),  LiNbO.3  (532  nm),  KTP  (532  nm),  and  tellurite  glass 
(532  nm).  For  materials  where  two  values  of  ft  were  obtained, 
an  average  ft  was  determined. 

We  now  use  the  newly  assigned  bandgap  energies  E'g  to 
recalculate  n2.  This  results  in  somewhat  better  agreement  with 
(3),  however,  there  is  still  a  sign  discrepancy  for  BBO  at 
266  nm,  and  the  experimental  n->  at  266  nm  for  calcite  is  23 
times  larger  than  predicted.  If  the  bandgap  for  BBO  is  further 
increased  from  6.45  eV  to  6.8  eV  (a  5%  increase),  the  n2 
datum  at  266  nm  can  be  made  to  agree  with  (3)  and  (4).  This 
only  leads  to  a  small  discrepancy  with  the  ft  values  at  266  nm 
(ft  —  0.77  cm/GW  rather  than  the  experimentally  observed  0.9 
cm/GW)  and  at  355  nm  (  ft  —  0.066  cm/GW  rather  than  the 
experimentally  observed  0.01  cm/GW).  This  shift  eliminates 
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Fig.  6.  Spectrum  of  the  scaled  experimental  values  of  13.  The  solid  line  shows  the  prediction  of  the  two-parabolic  band  model. 


the  discrepancy  in  sign  for  n2  at  266  nm  but  the  resulting 
value  for  n2  at  355  nm  does  not  show  an  enhancement  around 
the  new  2PA  edge  (which  now  lies  near  355  nm). 

The  discrepancies  for  CaCOa  are  even  more  significant. 
Shifting  Eg  to  E'g  ~  6.8  eV  (a  15%  increase)  gives  excellent 
agreement  for  the  0  values  at  both  266  nm  and  355  nm, 
however,  the  experimental  value  of  n2  at  266  nm  is  23  times 
larger  than  predicted.  In  order  for  calcite  data  to  match  the 
predicted  dispersion,  a  value  of  E'g  ~  8  eV  is  needed  (a  further 
increase  of  16%).  This  would  also  place  the  355-nm  datum 
below  the  2PA  edge  where  0  should  vanish.  However,  the  very 
small  value  of  0  at  355  nm  of  0.018  cm/GW  is  consistent  with 
indirect  gap  2PA.  For  this  value  of  0  to  be  due  to  direct  gap 
2PA  would  require  355  nm  to  lie  very  close  to  the  2PA  edge 
(note  that  this  is  also  true  for  BBO  at  355  nm  where  the  small 
experimental  value  of  0  is  consistent  with  indirect  gap  2PA). 
If  we  assume  that  calcite  is  an  indirect  gap  material  and  that 
E'  ~  8  eV,  the  ri2  data  closely  follow'  the  predicted  scaling 
relations  and  the  value  of  0  at  266  nm  is  predicted  to  be  0.33 
cm/GW  as  opposed  to  the  experimental  value  of  0.8  cm/GW, 
i.e.,  within  expectations. 

Fig.  7  shows  the  results  of  rescaling  the  data  using  E'g, 
plotted  with  the  predicted  dispersion  curve,  Cz ■  The  values 
used  for  Eg  are  listed  in  Table  I  (the  values  for  BBO  of  6.8 
eV  and  for  calcite  of  8  eV  were  used).  The  data  of  Fig.  3  are 
shifted  up  and  to  the  left  by  this  increase  in  bandgap  energy. 


We  now  obtain  better  agreement  for  n2  with  the  two-parabolic 
band  calculation  around  the  2PA  edge  where  some  data  show 
the  predicted  enhancement,  whereas  in  Fig.  3,  the  scaled  data 
remain  relatively  constant  throughout  the  measured  spectrum. 
The  largest  disparity  is  for  BBO  at  355  nm  which  is  a  factor 
of  4.5  lower  than  the  theory.  Most  data  falls  within  a  factor 
of  two  of  the  prediction.  Fig.  8  shows  a  log-log  plot  for  the 
rescaled  data  where  we  now  find  much  of  the  data  straddling 
either  side  of  the  predicted  E line  (all  data  now  lie  on  the 
graph). 

Rescaling  the  bandgap  energy  in  accordance  with  2PA  mea¬ 
surements  has  precedence.  Frohlich  and  Staginnus  [25]  have 
previously  shown  bandgap  reassignment  in  alkali  bromides  by 
using  two-photon  spectroscopy.  Bosio  and  Czaja  [26]  found 
that  using  a  bandgap  energy  of  9.1  eV  fits  their  analysis  of 
the  Urbach  edge  [27]  in  amorphous  silicon  dioxide  (aSi02). 
This  is  nearly  identical  to  our  Eg  of  9.0  eV  found  by  shifting 
the  2PA  data.  Ching  and  Xu  [28]  calculated  the  band  stmeture 
of  KTP  and  obtained  a  direct  gap  of  ~4.9  eV,  which  is  close 
to  our  E'g  of  ~4.6  eV.  Sapphire,  although  an  indirect  gap 
material,  may  be  approximated  as  a  direct  gap  insulator  as  the 
energy  difference  between  indirect  and  direct  gaps  is  small 
[29].  Measured  experimental  values  of  the  bandgap  in  sapphire 
range  from  8.7  eV  to  9.4  eV  [29],  Our  value  of  E'g  —  8.7  eV 
is  within  this  range.  Little  information  on  the  bandgap  of  BBO 
exists,  however,  BBO  and  LiB30.5  (lithium  triborate  or  LBO) 
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Fig.  7.  Plot  of  the  dispersion  of  n  rescaled  with  l'[t .  as  a  function  of  hujj  E'g . 


have  similiar  linear  optical  properties,  i.e.,  refractive  index 
and  transmittance  window.  The  electronic  structure  of  LBO 
was  calculated  by  Xu  and  Ching  [30]  from  which  a  direct  gap 
of  7.4  eV  was  found.  Because  of  the  structural  and  optical 
similiarities,  our  value  of  E’  ~  6.8  eV  for  BBO  appears  to 
be  reasonable. 

IV.  Conclusion 

We  have  extended  measurements  of  the  bound-electronic 
nonlinear  refractive  index,  n2.  and  two-photon  absorption  co¬ 
efficient,  ft,  to  wide  bandgap  dielectric  materials  to  frequencies 
available  from  the  harmonics  of  the  Nd:YAG  laser,  i.e.,  1064, 
532,  355,  and  266  nm.  The  purpose  of  this  study  was  both 
to  extend  the  database  and  to  determine  the  validity  of  the 
very  successful  two-parabolic  band  model  that  predicts  these 
nonlinearities  in  semiconductors  [6],  [7],  [9],  [19].  That  theory, 
which  predicts  the  nondegenerate  (and  degenerate)  nonlinear 
absorption  spectrum  and  used  a  Kramers  Kronig  analysis  to 
obtain  the  dispersion  of  n->.  gives  the  bandgap  scaling  and 
frequency  dependence  of  n2  and  ft  correctly  for  a  wide  variety 
of  materials.  It  also  gives  values  of  n2  for  wide  bandgap 
materials  to  within  a  factor  of  2  to  3  for  frequencies  far  below 
the  band  edge.  This  study  extends  the  frequency  range  of  data 
much  nearer  to  the  band  edge  in  these  materials.  One  of  the 
difficulties  in  determining  the  validity  of  the  two-band  model 
is  that  it  is  based  on  a  single  direct  bandgap.  In  many  of  the 


dielectric  materials  studied,  this  is  a  poor  approximation  for 
the  linear  behavior.  What  we  determine  is  that  if  the  bandgap 
is  defined  as  the  linear  absorption  cutoff  energy,  agreement  is 
relatively  poor,  and  even  the  sign  of  no  is  incorrectly  predicted 
for  two  of  the  materials  studied  at  266  nm.  However,  another 
method  previously  used  to  define  the  bandgap  is  the  energy 
of  the  two-photon  absorption  edge.  By  forcing  our  measured 
values  of  ft  to  fit  the  two-parabolic  band  model  by  allowing 
Eg  to  increase,  better  agreement  is  obtained  for  no  although 
significant  discrepancies  remain  for  BBO  and  calcite.  A  further 
increase  in  bandgap  energy  for  BBO  and  calcite  to  reach  a 
compromise  agreement  with  theory  for  the  n2  and  ft  data 
reduces  these  discrepancies.  These  new  bandgaps  ( E'g )  can 
be  considered  effective  direct  gaps  for  each  material.  The 
bandgap  energy  values  obtained  are  found  to  be  consistent  with 
previous  studies.  We  also  note  that  this  procedure  is  incomplete 
in  that  we  only  have  2PA  data  for  those  materials  with  bandgap 
energy  low  enough  to  exhibit  2PA  at  266  nm. 

This  two-parabolic  band  model  is  clearly  an  oversimplifi¬ 
cation  and  more  sophisticated  models  are  needed  to  give  the 
details  of  the  dispersion  of  no  and  spectrum  of  ft.  For  example, 
[31]  showed  the  existence  of  an  anisotropic  2PA  coefficient 
can  be  explained  by  including  the  effects  of  higher  bands, 
and  measurements  showed  good  agreement  [32],  Clearly, 
anisotropy  in  ft  also  leads  to  anisotropy  in  no.  However,  the 
overall  agreement  within  factors  of  2  or  3  for  such  a  large 
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class  of  materials  over  such  a  large  frequency  range  (1064 
nm  to  266  nm)  for  a  greater  than  102  change  in  magnitude  of 
rt2  shows  that  the  basic  physics  of  the  interactions  are  being 
included  here,  and  that  the  details  of  the  band  stucture,  while 
playing  a  significant  role,  do  not  play  the  dominant  role  in 
determining  these  nonlinearities. 
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